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The data of the experiment of Schiller et al., Phys. Rev. Lett. 77(1996) 2933, are 
alternatively evaluated using the maximum likelihood estimation. The given data 
are fitted better than by the standard deterministic approach. Nevertheless, the 
data are fitted equally well by a whole family of states. Standard deterministic 
predictions correspond approximately to the envelope of these maximum likelihood 
solutions. 



1 Introduction 

Quantum state provides the complete information about quantum systems. Recently 
quantum tomography has been devised for prediction of quantum state on the basis 
of homodyne detection with rotated basis of quadrature operators [1,2,3]. The tech- 
nique has been applied to analysis of realistic measurement and now, quantum state 
reconstruction is routinely used in various applications [4] . 

Nevertheless, potential problems of deterministic schemes has been reported. The 
positive definiteness of the reconstructed density matrix is not guaranteed within de- 
terministic data inversion yielding some nonphysical predictions. Positive definiteness 
can be preserved using information theory [5]. The approach based on the maximum 
likelihood (MaxLik) estimation is closely related to the standard treatment. Instead of 
the question: "What quantum state is determined by these data?" the question consis- 
tent with quantum theory reads: "What quantum state seems to be most likely ?" The 
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general formulation of the MaxLik problem was given in [6,7]. The extremum density 
matrix is given by nonlinear operator equation 

Rip)p = p, (1) 

where 

R = Y,—^i, Pii = Mpfl^). (2) 
. Pii 

Here Ilj represents the projectors corresponding to (in general unsharp) nonorthogonal 

measurement. The mcasiircd relative frcqiicncics arc denoted here as /; = 1. The 

state may be reconstructed on the subspace where the projectors provide the resolution 
of identity operator .R = 1. Although the form of equation (1) suggests an iterative 
solution, it is not easy since the iterations need not be convergent in general. Solution 
simplifies significantly provided that the projectors Hi commute. The diagonal elements 
of density matrix in common commuting basis may be estimated very effectively, as 
showed by Banaszek [8,9]. This numerical approach will be used here for evaluation of 
the diagonal elements of density matrix for the experiment reported by Schiller et. al 
[3]. All the necessary steps of general reconstruction scheme will be demonstrated on 
this example. 



2 Reconstruction 

The diagonal elements of density matrix only will be reconstructed. Data corresponding 
to the random-phase homodyne detection [10] are sufficient for this purpose. Projectors 
enumerated by position are given as 

ti{x) = ^j\e\x,0){x,e\. (3) 

They commute and are complete on the full interval x G (—00,00). Nevertheless, any 
realistic measurement will register only a finite sampling of discrete decomposition. 
Denoting the position of a particular bin in x-coordinate as Xi, the projectors in number 
state basis read 

00 n / \ 

n(xO =AxJ2J2 - VT-' (I) \n){n\, (4) 

Here r] denotes the efficiency of counting of photoelectrons, enumerated by index k 
and Aa; is the width of the bin. This detection of discretized quadrature components 
reproduces the identity operator as 

^ = ^%0- (6) 



Reconstruction of diagonal elements. 



3 



Diagonal elements of R in number state basis shows the subspace where this operator 
equals approximately to identity. Here the reconstruction may be done applying the 
MaxLik procedure. Relative entropy (log of likelihood function) 



shows how the given state approaches the ideal condition pu = fi. Provided that this 
condition is met, the relative entropy equals to zero. In the following the relative 
entropy will be evaluated for various states and given data fi. For comparison, the 
relative entropy will be expressed in % with respect to the entropy of measured data. 



The data corresponding to measurement of coherent and squeezed state will be con- 
sidered explicitly. All the calculations has been done for an overall efficiency rj = 0.85 
here. The Fig. 1 plots the data measured by random-phase homodyning and recon- 
structed data using deterministic and MaxLik approaches. The subspace for reconstruc- 
tion is deliminated by the dimension about Uedge = 50 as follows from the plot of the 
diagonal elements R{n) in the Fig. 2. The reconstruction of diagonal elements of den- 
sity matrix in number state basis is plotted in the Fig. 3. Deterministic reconstruction 
yields slightly nonphysical results indicated by several negative diagonal elements ob- 
tained. A typical MaxLik estimation of diagonal elements is plotted in the upper right 
panel. However, the estimates depend on the starting point of the iteration process. 
The procedure has been repeated 100 times with different starting points. The average 
of all estimates is plotted in left lower panel of the Fig. 3. The MaxLik estimations 
fits the measured data obviously better. While the relative entropy for the determin- 
istic "state" is about K{p/ f) w 1.6% of the value S = 4.717, the relative entropy of 
the MaxLik estimates fluctuates around the value 0.31% of the entropy S only. The 
histogram of relative entropies is plotted in the left upper panel of the Fig. 4. Sig- 
nificantly, all these states fit well the input data and there is no observable difference 
in fitted data statistics. Other panels show the uncertainty of moments for various 
MaxLik estimations. Moments are normalized with respect to the average of all the 
MaxLik estimates (denoted by index av)- The x-coordinate represents the deviation 
of the moments {h'^) / {n'^) av — 1 for /c = 1, 10, 50. All the low moments are estimated 
very sharply within the accuracy 10~^%. However, the 50-th moments already fluctuate 
within 30%. The average number of particles is estimated here as fiMaxLik ~ 29.6. The 
value obtained from the deterministic approach is n^et ~ 25.6, i.e. 30.1 considering 
efficiency r}. 

Similar analysis may be done for the data corresponding to squeezed state. The 
histogram of random-phase homodyne detection and reconstructed data are plotted in 
the Fig. 5. The dimension of subspace for reconstruction is about nedge = 100 (not 
plotted here). The reconstructed diagonal elements are shown in the Fig. 6. The 
deterministic reconstruction has been considered on too small subspace as seen on the 
upper left panel. The MaxLik estimation has been done 100 times on a 75 dimensional 
subspace, becaiise the higher dimension is out of the range of our program. While 
deterministic estimation is characterized by the relative entropy Kgg{p/ f) 15% of 




(7) 



'?(/) = - Ei/i In /i- 
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the value S = 4.02, the relative entropy of MaxLik estimates fluctuating around the 
value 0.15% are considerably better. The histogram of relative entropies is plotted in 
the upper left panel of the Fig. 7, the other panels show the uncertainty in moments 
of various MaxLik estimates. The value of moments is again related to the averaged 
MaxLik estimate as in the Fig. 4. 

3 Conclusion 

The MaxLik procedure fits the measured data better than the deterministic scheme, but 

the results are not single valued. Instead of a single state predicted by the deterministic 
scheme, there is a family of states fitting the data equally well. Hence the MaxLik 
state reconstruction is more uncertain in comparison to the deterministic prediction. 
However, this uncertainty corresponds to the probabilistic nature of quantum theory. 
Measured data do not determine an unique state and the reconstruction has to take it 
into account. 
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Fig. 1: Homodyne detection and reconstructed data for coherent state 



Fig. 2: Diagonal elements of decomposition of identity 



Fig. 3: Reconstruction of diagonal elements for coherent state 



Fig. 4: Histogram of relative entropies and moments of particle number operator for various 
MaxLik estimates. 



Fig. 5: Homodyne detection and reconstructed data for squeezed state 



Fig. 6: Reconstruction of diagonal elements for squeezed state 



Fig. 7: Histogram of relative entropies and moments of particle number operator for various 
MaxLik estimates 
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